Symplectic tomography of ultracold gases in tight-waveguides 
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The phase space is the natural ground to smoothly extrapolate between local and non-local cor- 
relation functions. With this objective, we introduce the symplectic tomography of many-body 
quantum gases in tight-waveguides, and concentrate on the reduced single-particle symplectic to- 
mogram (RSPST) whose marginals are the density profile and momentum distribution. We present 
an operational approach to measure the RSPST from the time evolution of the density profile after 
shutting off the interactions in a variety of relevant situations: free expansion, fall under gravity, 
and oscillations in a harmonic trap. From the RSPST, the one-body density matrix of the trapped 
state can be reconstructed. 
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At low densities, ultracold bosonic gases exhibits 
universality. The interatomic interactions are then 
well described by the Fermi-Huang pseudo-potential, 
parametrized by a the 3D s-wave scattering length a s . 
If such gases are further confined in tight- waveguides, 
whenever the transverse excitation quantum huj± is 
larger than the longitudinal zero point and thermal ener- 
gies, the system effectively becomes one-dimensional [![. 
The interparticle pseudo-potential is then a simple delta 
function, so the system is well approximated by the Lieb- 
Liniger model 0. Moreover, the strength of the interac- 
tion as a function of a s exhibits a confinement-induced 
ID Feshbach resonance (CIR) allowing to tune the 

ID coupling constant g\o from — oo to +oo and to reach 
both weak and strongly interacting regimes [1] . As a con- 
sequence, paradigmatic examples of the Bose-Fermi du- 
ality have been explored such as the Tonks- Girardeau gas 
H, in which the strongly repulsive interactions between 
bosons leads to an effective Pauli exclusion principle @. 
In this regime, the system undergoes fermionization, all 
local correlation functions being identical to those of the 
spin-polarized ideal Fermi gas. Actually, the Fermi-Bosc 
duality comes into play even with finite interactions Q- 
However, quantum statistics invariably imposes an un- 
derlying signature manifested when looking at non-local 
correlations such as the momentum distribution or the 
one-body density matrix | 1 1 [EI1I3. A nat- 
ural ground to smoothly extrapolate between local and 
non-local correlations is the phase space. In this paper, 
we shall undertake the description of ultracold gases in 
tight-waveguides by means of the quantum tomographic 
technique [14(. We show that after shutting off the in- 
teractions in the system, the time evolution of the den- 
sity profile governed by a quadratic Hamiltonian is tanta- 
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mount to the knowledge of the (reduced) symplectic to- 
mogram, from which the initial one-body density matrix 
of the trapped state can be reconstructed. This includes 
relevant experimental situations for ultracold atoms in 
waveguides such as free expansion, dynamics falling un- 
der gravity, and time-evolution in a harmonic trap. 

Symplectic tomography. In the symplectic tomogra- 
phy probability representation, first introduced in [14j . 
the wave function ^(z) or the density matrix p(z, z') 
can be mapped onto the standard positive distribution 
W(X, [i, v) of the random variable X depending on two 
real extra parameters, \x and v. The map is given by the 
formula 

W(X,n,v) = —Jr-r / p(v')^"^ (z - z,) <1^, 
ZttTi\v\ J 

(1) 

which is the fractional Fourier transform [151 ] of the den- 
sity matrix. The map is invertible so that the density 
matrix can be expressed in terms of the tomographic 
probability representation as follows, 

PM = ^. J W(X,^^)e< X -^)dXdp. 

(2) 

The expression in Eq. |T]) admits an affine invariant form 

W{X,n,v) = Tt\p6{X-f£-vp)], (3) 

where the density operator is denoted by p, and z, p 
are the operators of position and the conjugate momen- 
tum respectively. From Eq. ([3]) some properties of 
the tomogram W(X, fi, v) are easily extracted. First, 
the tomogram is a normalized probability distribution, 
f W(X, fi, u)AX = 1, if the density operator is accord- 
ingly normalized (i.e. Trp — 1). Moreover, the tomogram 
satisfies the homogeneity property [l8[ 

W(XX, Xfx, Xv) = ±rW(X, (i, u), (4) 
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inherited from that of the delta function in the definition, 
Eq. ([3]). This equation provides a link with optical to- 
mography (p = cos 9, v = sin 9) and Fresnel tomog- 
raphy ([a — 1) [20(. We further notice that the tomogram 
can be related to the Wigner function W(z,p) [2l| . 



W(X,n,v)= W(z,p)5(X-ftz-vp) 



dzdp 
2irh 



(5) 



as its Radon transform, providing a clear interpreta- 
tion of the tomogram W(X, fi, v). One has the line 
X = pz + vp in phase space, which is given by equating 
to zero the delta- function argument in Eq. ©. Alterna- 
tively, the parameters p and v can be expressed in the 
form s cos 6, v = s~ 1 sin 9. Here s > is a real squeezing 
parameter and 9 is a rotation angle. Then the variable 
X is identical to the position measured in the new refer- 
ence frame in the phase-space. The new reference frame 
has new scaled axis sz and s _1 p and after the scaling 
the axis are rotated by an angle 9. Thus the tomogram 
implies the probability distribution of the random posi- 
tion X measured in the new (scaled and rotated) refer- 
ence frame in the phase-space. The remarkable property 
of the tomographic probability distribution is that it is 
a fair positive probability distribution and it contains a 
complete information of the state. Indeed, the density 
operator p can be expressed in terms of the tomogram as 



i(X — fiz — up) 



dXdpdv. (6) 



The tomographic map can be used not only for the de- 
scription of the state in terms of a probability distribu- 
tions, but also to describe its evolution (quantum transi- 
tions) by means of the standard real positive transition 
probabilities (alternative to the complex transition prob- 
ability amplitude), i.e., in tunnelling dynamics |23j . 

Many-body ID gases. We next focus on effectively one- 
dimensional many-body systems, described generally by 
a wavefunction ij)(z\, . . . , Zjsr) or alternatively the N-body 
density matrix Pn(z±, . . . ,zn, z[, . . . , z' N ). Introducing 
the notation z = z±, . . . , zn, and similarly for {X , u, fi}, 
one finds the many-body tomogram carrying out the 2N- 
dimensional integral, 



1 



p N (z,z') 



E 3 



xe 



dzdz 



(7) 



Let us define the reduced single-particle symplectic to- 
mogram (RSPST) as 



W(X 



N 

,H,v)= / Y\dX 3 W N {X,^u). (8) 

J 3=2 



From the basis invariance property of the trace, it 
follows that alternatively one can find W(X, p, v) 



through Eq. ([lj, using the reduced single- 
particle density matrix (RSPDM) p(z, z') = 
J dz2 ■ ■ ■ dzNip{z, Z2, ■ ■ ■ , zn)iP(z', Z2, ■ ■ ■ , zn)* • Note that 
we choose the normalization J dzp(z, z) = 1. Further 
notice the normalization condition J W(X, p, v)dX = 1, 
and that the symplectic tomogram satisfies the following 
two marginals of the Wigner function, 



W(«,l,0) = n(z), 



W(p,o,i) = e (p) ) 



(9) 



where n{z) — p(z, z) is the density profile, and g(p) = 
(2irh)~ 1 J dzdz' p(z, z')e lp ( z ~ z ^ h the momentum distri- 
bution. Therefore, the reduced tomogram extrapolates 
smoothly between the (local) density profile and the 
(non-local) momentum distribution performing a rota- 
tion in phase space parametrized by (p, v). 

Measuring the RSPST. A direct experimental measure- 
ment of the matter- wave symplectic tomograms has been 
eluded so far. To close this gap, we shall consider dif- 
ferent situations to implement physically the symplectic 
tomography representation. We shall next focus on the 
Lieb-Liniger model [2J which accurately describes ultra- 
cold atom vapors strongly confined in waveguides [J H, 0] . 
The Hamiltonian is that of N trapped bosons with pair- 
wise delta interactions, 



N 



H 



LL 



(10) 



l<i<j<N 



where Hi = 

1 Ira az 



+ V(zi) is the single-particle Hamil- 
tonian, and V(z) denotes the trapping potential. The 
value of g\D is a function of a« with a confinement- 
induced ID Feshbach resonance [l|, which allows to 
tune g\D from — oo to +oo by shifting the position of a 
3D Feshbach scattering resonance via an external mag- 
netic field 01. 

In what follows we shall consider the time-dependence 
of the coupling strength to be gmit) = t) (where 

<d(t) is the Heaviside step function), so that the interac- 
tions are turned off for t > 0. The RSPST can then be 
related to the dynamics governed by an external poten- 
tial. 

If for t > the trap is also turned off the Hamilto- 
nian becomes the kinetic energy operator, and the sys- 
tem undergoes free expansion. The free time evolution is 
governed by the kernel (2f| 

fC (z,z',t\z ,z' ,0) = (z' \tf (t)\z')(z\U Q (t)\zo) (11) 

where the matrix elements (propagator) of the free evo- 
lution operator Uo(t) read 



(z\U (t)\zo) 



rn 



2niht 



(12) 



so that 



p (z,z',t)= dz a dz' K. (z,z',t\z ,z' ,0)p(z a ,z' ). (13) 
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It follows that the RSPST can be obtained from the free 
propagation of the density profile no(z, t) = po(z, z, t), 



\ p, P 



(14) 



The expansion takes place here in one dimension rather 
than three, this is, while keeping the radial confinement 
on. We further notice that if interactions are kept, im- 
portant deviations from the ballistic expansion may occur 
[26L [27J , and the Lieb-Liniger dynamics blurs the infor- 
mation of the initial state [28|, [29| • 

A careful analysis in the presence of the gravitational 
potential V(z) — mgz, whose propagator can be related 
to the free one in Eq. (fT2l as (z\U~o(t)\zo}exp(— img(z + 
z )t/2h-im 2 g 2 t 3 /24h)lgl^, allows to obtain the RSPST 
in an analogously way from the density profile, with the 
replacement z — > X/p — gm 2 v 2 /2p 2 . Therefore, the ex- 
pansion both in a horizontal or vertical waveguide allows 
to reconstruct the reduced tomogram and hence the one- 
body density matrix of the initial state. 

Alternatively, the RSPST can be measured in a har- 
monic trap without releasing the gas, just by shutting off 
the interactions. In this case the tomographic kernel in 
Eq. ([I]) is implemented using the propagator [251 ] 



(z\Utrap{t)\z ) 



2irihsmujt 



g z 51 \ z + z al l n Biiiui , 



whence it follows that 



W(X,p,v) = -n t 



1 / X 1 i rruvv 



trap I » j 



A U) 



tan 



(15) 



(16) 



with A = \J p 2 + m 2 uj 2 v 2 



At variance with the oper- 
ation approach based on free expansion, the dynamics 
is periodic, and it is possible to reconstruct the RSPST 
from the density profile along one cycle, T = 2it/uj. 

Reconstruction of the density matrix. The reconstruc- 
tion of a quantum state from the dynamics of wavepack- 
ets in a potentail V(z) has been successfully addressed by 
Leonhardt and coworkers within the optical tomography 
[13, [H|. Moreover, the Wigner function has been exper- 
imentally measured by looking at the time-evolution of 
the density profile of Helium atoms [32l |: and recently, 
a scheme for neutron wavepacket tomography has been 
proposed [11]. 

In the following, we provide the explicit expressions 
for the density matrix corresponding to the situations 
discussed in the previous section. Note that using 
Eqs. (O and (| 14[) . the free time evolution of the 
density profile can be written in the invariant form, 
no CM) = (z\U (t)p(%(t)\z) = Tr [p5(z - z - tp/m)], 
which suggests that this type of reconstruction can be 
applied to other observables different from p. More- 
over, from Eq. ^ and (|14p . defining X = kz, using 
the Baker-Campbell-Hausdorff formula and the Jacobian 



\d(X, v, /i)/<9(k, z,t)\, the inverse reads 



2nm 



n () 



{z,t)\K\e l < z - i -^ p )dKdzdt. (17) 



Let as denote the Fourier transform of the density profile 
by ho(k,t) = -j== f n (z,t)e~ lkz dz. In the wavevector 
^-representation one finds that the density matrix of the 
initial state is given by 



p(k,k') = 



27T7TJ 



no (k — k , t) | k — k'\e 



h(k 2 -k' 2 )t 



: dt,(18) 



which is analogous to the equation derived in [3l| for 
wavepacket state reconstruction. The Riemann-Lebesgue 
lemma limits in practice the range of the integral. Al- 
ternatively, we note that the RSPDM can be diagonal- 
ized [13] as p(z,z';t = 0) = Ylj ^j<Pj{ z )<Pj i z ') m terms 
of the orthonormal natural orbitals fj(z) with occupa- 
tion numbers Xj > satisfying j p{z, z')ipj(z')dz' = 

Xjipj(z) and J2jLi -\? = 1- Under free evolution, 
having set up gm{t > 0) = 0, the density pro- 
file reads no(-M) = Xj\ipj(z, t)\ 2 with tf>j(z,t) = 

(271-)- 1 / 2 / dk(f j (k)e ikz - lhk2t / 2m , which using Eq. ([15]). 
leads to the density matrix p(k 7 k') = Y]j Xj(pj{k)(p*(k') 

When the density profile corresponds to the expan- 
sion falling under gravity, the above expression Eq. 
(JTSJ holds with the definition ?1o(fc,i) = ^- J n (z — 
gt 2 /2,t)e- lkz dz. 

The dynamics in a harmonic trap after shutting off the 
interactions similarly allows to reconstruct the density 
matrix operator according to 



P = 



2nm 



J n trap (z,i)| K |e 4K ( z - cos " ti -^ i5 )dKdzdt, 



(19) 



or in the wavenumber representation, 
ft 



p(k,k') 



2wm J cos 2 cot 



1 „ fk-k' 



COS LUt 



x\k-k'\e i& ^( k2 - k ' 2 Ut, 



(20) 



where the caustic at half period is a well-known property 
of the Radon transform. We note also that the obtained 
results can be put in the compact form 



P : 



27T77J 



(z,t)\K\e^ z - z ^d K dzdt, (21) 



where z(t) = U h zUh, and nn(z,t) is the density profile 
evolving under the action of Uh, the evolution operator 
corresponding to a quadratic Hamiltonian H of free mo- 
tion, gravitational fall, and harmonic oscillation. This 
equation suggests that the reconstruction of the density 
matrix can be generalized for time dependent single par- 
ticle quadratic Hamiltonians. 
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Conclusions. In conclusion, we have introduced the 
reduced symplectic tomography of many-body quantum 
gases confined in tight-waveguides. By performing a ro- 
tation in phase space, a smooth extrapolation is possible 
between the (local) density profile and the (non-local) 
momentum distribution. Moreover, a simple procedure 
has been given to experimentally measure the tomogram, 
namely, by suddenly shutting off the interactions and 
registering the time evolution of the density profile in 
a variety of situations: free expansion, expansion under 
gravity, and periodic motion in a harmonic trap. We note 
that such dynamics can only be implemented if the in- 
terparticle interactions are negligible. Else, the so-called 
cusp condition imposed on the wavefunction by the short- 
range pseudo-potential affects the density profile [13, H§| 
blurring the information of the initial state. In addition, 
it should be clear that even though the density profile of 
duals systems within the trap is exactly the same, once 
the interactions are switched off, the time evolution will 
differ provided that the momentum distribution is un- 
like in each system and therefore our method account 



for a proper description of different quantum statistics. 
The one-body density matrix of the state in the trap can 
then be reconstructed. As long as the dynamics is gov- 
erned by a quadratic Hamiltonian, the reconstruction is 
expected to be possible even for time-dependent poten- 
tials. We close by noting that higher order correlations of 
the trapped gas can be inferred from the time evolution 
of the density profile [3(| . 
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